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Abstract — Binary LDPC coded relay systems have been well 
studied previously with the assumption of infinite codeword 
length. In this paper, we deal with non-binary LDPC codes 
which can outperform their binary counterpart especially for 
practical codeword length. We utilize non-binary LDPC codes 
and recently invented non-binary coding techniques known as 
multiplicative repetition to design the low-rate coding strategy for 
the decode-and-forward half-duplex relay channel. We claim that 
the proposed strategy is simple since the destination and the relay 
can decode with almost the same computational complexity by 
sharing the same structure of decoder. Numerical experiments 
are carried out to show that the performances obtained by non- 
binary LDPC coded relay systems surpass the capacity of direct 
transmission and also approach within less than 1.5 dB from the 
achievable rate of the relay channels. 

Index Terms — non-binary low-density parity-check code, mul- 
tiplicative repetition, low-rate code, rate-compatible code, decode- 
and-forward, relay channel 

I. Introduction 

Over recent years, cooperative communication has been 
studied extensively in order to increase the capacity of wireless 
networks H], 121, 131. In cooperative networks, terminals 
act as relays which can forward the received information to 
destinations or terminals. It is well known that spatial diversity 
is the promising technique for improving the quality and 
reliability of wireless links IH. The spatial diversity can be 
obtained by adopting multiple antennas in a transmitter or/and 
a receiver 15), Q. Relays in a cooperative network form a 
virtual multiple antenna system and hence allow a single- 
antenna user to transmit the information through the different 
and independent paths. Without employing multiple antennas, 
the spatial diversity still can be achieved via relaying. There- 
fore, the concept of relaying is applicable to mobile terminals 
which cannot support multiple antenna. However, cooperative 
communications with multiple antennas also have been studied 
to achieve higher capacity Q, IS) but this is beyond the scope 
of this paper. 

A single relay channel is an elementary component of 
a cooperative network The single relay channel simply 
consists of three terminals which are a source, a relay and a 
destination. It is theoretically shown that the achievable rate 
of a single relay channel is much higher than the capacity 
of direct transmission ITOl . Many protocols, such as amplify- 
and-forward ITTI or compress-and-forward lT2l . have been 
developed for processing the signals in relay channels. Previ- 
ous works have shown that the decode-and-forward protocol 
achieves better performance compared with other protocols 
when a quality of source-relay link is good lT3l . lT4l . By using 



the decode-and-forward protocol, a relay is capable to decode 
the received signal and forward the decoded information to 
the destination. 

At the present moment, half-duplex relaying is regarded 
more practical since full-duplex relaying cannot be imple- 
mented efficiently. The operation in a half-duplex scheme is 
divided into two modes named as broadcast (BC) and multiple 
access (MAC) modes. In BC mode, the source transmits signal 
to both the relay and the destination. The relay does not have 
its own information to transmit. Thus the relay only listens 
to the source in BC mode. In MAC mode, both the relay 
and the source transmit their signals to the destination. In this 
paper, we focus only on the half-duplex relay channel with 
the decode-and-forward protocol and restrict ourselves to the 
case of a single relay. 

Several binary LDPC (BLDPC) coding strategies have 
been developed to enhance the performance of a decode- 
and-forward half-duplex relay channel ITSl . lfT6l . The main 
challenging problem in this field is that one needs to design 
a coding strategy that works at both source-relay and source- 
destination links. In other words, a code-structure is necessary 
to exhibit excellent performance at two different signal to 
noise ratios. After designing the code properly, the code is 
distributed to the source, the relay and the destination. We 
briefly introduce some prominent works in this research field. 
Chakrabarti et al. optimized the degree profile of BLDPC 
codes for half-duplex relay channel by mean of density evo- 
lution ITtI . Razaghi and Yu developed a technique called bi- 
layer BLDPC codes QH in which overall graph of an BLDPC 
code is subdivided into lower and upper parts corresponding 
to codes for source-relay and source-destination links. In lT6l . 
Li et al. proposed the distributed BLDPC code using a rate- 
compatible structure and also developed a method which can 
accurately predict the performance of BLDPC codes for a relay 
channel. Recently, Cances et al. extended the idea of bi-layer 
and designed the BLDPC codes that work closely to theoretical 
limit lT9l . We note that all works mentioned above focused 
on BLDPC codes. 

Although the significant progress has been made to design 
the distributed LDPC codes for the half-duplex relay channel. 
However, there still exist some problems that can be addressed 
as follows; 1) The design of optimized LDPC codes for the 
relay channel requires a task to obtain such an optimized de- 
gree profile. 2) The sophisticated degree profile of LDPC code 
obtained from an optimization process makes the hardware 
implementation quite complicated. 3) LDPC codes with the 
optimized degree profile require very large codeword length 



2 



(e.g. 10^ bits) to perform close to the theoretical limit. The 
large codeword length leads to the transmission latency and 
complex decoders, which are not preferred in the real world 
communication. Moreover, the good performance of optimized 
LDPC code are not guaranteed when the codeword length is 
short or moderate. In this paper, we overcome these problems 
by considering non-binary LDPC (NBLDPC) codes. 

An NBLDPC code is defined by a sparse parity-check 
matrix defined over GF(2™), m > 1 ||20ll . Over point-to-point 
channels, it was shown in 1201, ||2T] that the NBLDPC codes 
can outperform their binary counterparts especially for the 
short and moderate codeword lengths. Therefore, the NBLDPC 
codes can be employed to improve the performance of prac- 
tical wireless communications since the short and moderate 
codeword lengths (e.g. a few thousand bits) are now currently 
used ED, El. Good NBLDPC codes defined over higher 
order Galois fields tend to have the regular degree profile 
(row and column weights of parity-check matrix are constant) 
l24l . The regular NBLDPC codes defined over GF(2™) with 
column weight 2 are empirically known [|24j as the best 
performing codes for 2™ > 64. Therefore, the first and second 
problems addressed in the previous paragraph can be alleviated 
since we can use NBLDPC codes with regular degree profile 
for transmission and hence do not need to optimize the degree 
profile of NBLDPC codes. 

At this moment, little progress has been made for an 
application of NBLDPC codes on relay channels but the afore- 
mentioned reasons mentioned earlier motivate us to design 
the distributed NBLDPC codes for relay channel. This paper, 
we develop the low-rate NBLDPC coding strategy for the 
decode-and-forward half-duplex relay channel. We apply the 
concept of multiplicative repetition ||251 . Il26l to design the 
distributed low-rate NBLDPC codes for the relay channel. 
Based on the proposed NBLDPC coding strategies, the relay 
and the destination can perform the decoding process with 
the same computational complexity even the coding rate of 
the destination is much lower than that of the relay. Finally, 
we show that the proposed strategy provides the relay systems 
with very good decoding performances especially for the short 
and moderate codeword lengths. 

The rest of this paper is organized as follows. Section II, 
we explain the half-duplex relay channel in time-division mode 
together with the decode-and-forward protocol. The NBLDPC 
codes are introduced in Section III. Section IV, we develop 
the simple NBLDPC coding strategy which applicable to low- 
rate relay channel. In Section V, we present the decoding 
performance of the proposed coding strategy. The conclusions 
are finally given in Section VI. 

II. Introduction to NBLDPC codes 

Before describing the proposed coding strategy, we first in- 
troduce NBLDPC codes, multiplicative repetition, puncturing 
and their decoding algorithm. 

A. NBLDPC codes 

NBLDPC codes are defined by sparse parity-check matrices 
defined over GF(2™), m > 1 ll20l . We can represent each 



element x E GF(2'") by using binary sequences of length 
TO bits. For transmitting over the binary input channels, each 
non-binary symbol in GF(2™) needs to be represented by a 
binary sequence of length m. For each to, we fix a Galois 
field GF(2'") with a primitive element a and its primitive 
polynomial tt. Once a primitive element a of GF(2™) is fixed, 
each symbol is given a TO-bit representation 1271 pp. 110]. 
For example, with a primitive element a £ GF(2'^) such that 
7r(a) = a"^ + a + 1 = 0, each symbol is represented as = 
(0,0,0), 1 = (1,0,0), a = (0,1,0), = (0,0,1), = 
(1,1,0), a"* = (0,1,1), = (1,1,1) and = (1,0,1). 
Assigning of binary sequence to each element depends on the 
primitive polynomial used to construct the field. Specifically, 
an NBLDPC code C over GF(2™) is defined by the null-space 
of a sparse M x N parity-check matrix H — {hi-) defined 
over GF(2'") 

C = {x e GF(2'")'^ I Hx = e GF(2")^^}, (1) 

where x = {xi, . . . , x^) is the codeword. The c-th parity- 
check equation for c ~ 1, . . . , M is written as 

/ici^i + hc2X2^ hc^XN = 0, (2) 

where h^, ■ ■ ■ , hc,^ S GF(2™) are the entries of c-th row 
of H. The parameter N is the codeword length in symbol. 
Throughout this paper (assuming H is of full rank), we denote 
the number of information symbols by K ~ N — M. Let n ~ 
mN and k = mK be the codeword length and information 
length in bit, respectively. The coding rate R of any NBLDPC 
code can be computed by i? = 1 — dy/d^ = K/N. 

An NBLDPC code is called (dv, c'c) -regular if the parity- 
check matrix of the code has constant column weight d^ 
and row weight dc- In this paper, only (c?v = 2, dc) -regular 
NBLDPC codes defined over GF(2^) are considered, since 
they are empirically known as the best performing codes. The 
dv,dc indicates the column and row weights of parity-check 
matrix, respectively. We can represent the parity-check matrix 
of an NBLDPC code by a Tanner graph with variable and 
check nodes |^8l. Each variable node and check node repre- 
sents a coded symbol and a parity-check equation, respectively. 
Therefore, the number of variable nodes and check nodes are 
equal to M and N, respectively. Figure [T] shows an example 
of the parity-check matrix of a (2, 3)-regular NBLDPC code 
defined over GF(2^) and the corresponding Tanner graph. 

Since, in this paper, we focus only on the binary transmis- 
sion, we describe how to transmit the NBLDPC coded symbols 
Xv by using the binary modulation scheme for v = 1,2, . . . , N 
and i = 1,2,..., to. At the v-th output of NBLDPC en- 
coder, each coded symbol Xy £ GF(2™) is mapped to the 
binary sequence of to bits {xy,i,Xv^2, ■ ■ ■ ,Xv,m) S GF(2)™ 
according to the primitive polynomial as described above. The 
obtained binary sequence {xv,i,Xy^2, ■ ■ ■ ,Xv,m) G GF(2)™ is 
then mapped to m modulated signals through the mapper for 
the transmission. 

B. Multiplicative Repetition 

More recently, Kasai et al. proposed an efficient method 
called multiplicative repetition to construct low-rate NBLDPC 
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Fig. 1. An example of (2, 3)-regular NBLDPC codes. The codeword length 
of this code is = 6. The number of parity symbol is M = 4. The upper 
part shows the parity-check matrix of size 4x6. The lower part shows the 
corresponding Tanner graph of this parity-check matrix but here we omitted 
edge labels for simplicity. Circle and square nodes represent variable and 
check nodes, respectively. 



codes ||251 . 1261 . Over point to point channels, the low-rate 
codes constructed by this method outperform the previously 
found low-rate codes. By using the NBLDPC code of rate 
Ri, we can easily construct code NBLDPC codes of lower 
rate R < Ri as described as follows. 

Let Ci denotes a NBLDPC code of length N and coding 
rate Ri. Since code Ci is used to constructed another codes, 
we refer to Ci as a mother code. A low-rate code C2 of 
length 2N and coding rate R2 = can be constructed 

as follows. We select N coefficients r^v+i, • • ■ , r2N randomly 
from GF(2'")\{0}. Note that we define = 1 e GF(2™) for 
V = 1, . . . ,N for simplicity of notation. We then multiplica- 
tively repeat the coded symbols from Ci with the coefficients 
to obtain the lower-rate code C2 as follows. 

C2 = {(Xi, . . . , X2n) e GF(2™)2^ I XN+v = rN+vXr,, 

for V ~ 1, . . . ,N, (xi, . . . , xn) e Ci }. 

In this way, we can construct a (2A^, K) code from an (N, K) 
code. 

The codes C3, C4, . . . , Ct of lower coding rates can also 
be constructed from code Ci through the multiplicative rep- 
etition process. We refer a parameter T as repetition pa- 
rameter. For T > 3, in a recursive fashion, N coefficients 
r(j-_i)7V-i-i: ■ • ■ 7 fTN are chosen randomly from GF(2™)\{0}. 
The code Ct of rate Rt = ^Ri can also be constructed as 
follows. 



CT = {(xi,...,XTiv)GGF(2'")™| 

X(_T-i)N+v = r^T-i)N+vXv, for V 
(xi, . . . ,X(T-i)jv) e Ct-i}- 



Therefore, we can easily construct {TN, K) code from {N, K) 
code in a recursive fashion. 

For example, we can construct C2, C3, . . . , Ct with coding 
rate 1/6, 1/9, . . . , 1/3T respectively when the mother code 
Ci is a (2, 3)-regular NBLDPC code. Figure |2] shows the 
block diagram of multiplicatively repeated encoder of Ct- We 
will show in the next section that the multiplicative repetition 



is applicable to design low-rate coding strategy for relay 
channels. 
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Fig. 2. Block diagram of multiplicatively repeated encoder of Ct- In this 
figure, P2S represents pai'allel to serial operation. 



C. Puncturing 

By puncturing a code, we obtain a higher-rate code. We 
cancel some parity symbols produced from the encoder of 
mother code C of rate i? = ^ to form the code C" of 
higher rate R' = where Np is the number of punctured 

symbols. We refer to the code C as the mother code. The 
puncturing plays an important role in constructing the rate- 
compatible codes. The rate-compatible codes are very useful 
in wireless communications which employ hybrid automatic 
repeat request (HARQ) protocol. The rate-compatible codes 
can flexibly adapt the coding rate according to the channel 
quality by adopting only one mother code. 

Instead of random puncturing, we find the location of 
variable node for puncturing by the concept called recoverable 
step 129] ■ This method helps the decoding process at R because 
the punctured variable nodes will receive nonzero messages 
at small number of iterations. We first search for the group 
of variable nodes with one-step-recoverable. After that, we 
then search for the group of variable nodes with two-step- 
recoverable, and so on. We therefore puncture variable nodes 
by the ascending order of recoverable steps. An example of 
puncturing on the (2, 3)-regular NBLDPC codes with one step 
recoverable is illustrated in Fig. [3] In this paper, the concept 
of puncturing by recoverable step will be used to obtain the 
rate-compatible NBLDPC codes for the relay channel. 




punctured 
node 



punctured 
node 



Fig. 3. An example of puncturing on (2, 3)-regular NBLDPC codes with one 
step recoverable. The codeword length of this code is Af = 6. The number 
of parity symbol is hi = 4. Circle and square nodes represent variable and 
check nodes respectively. We intend to puncture the variable node with one 
step recoverable first and then continue to puncture variable node with two 
step recoverable and so on. 
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D. Decoding algorithm for NBLDPC codes 

With some abuse of notation, we denote the w-th variable 
node by v. Let be the random variables with realization 
where represents the coded symbol and v = 1, 2, . . . , A^. 
Let be the random variables with realization y„ which is the 
received value from the channel Pr(F„|X„). The probability 
of transmitted symbol Pr(Xi,) is assumed to be uniform. We 
assume that the decoders at both relay and destination knows 
the channel transition probability 



Pr(X„=x|n=j/„),w==l,2, 



.,N, 



(3) 



for X G GF(2™). For memoryless binary-input output- 
symmetric (MBIOS) channels, the channel transition proba- 
bility is given as 

m 

Pr(X„ = x\Y^ = = nPr(X„,, = = (4) 

1=1 

where X^,i is the random variable of the transmitted bit, 
{xi,X2, ■ ■ ■ ^Xm) G GF(2)'" is the binary representation of 
X £ GF(2'") and the corresponding channel output and 
its random variable Ku.;. 

The decoding of NBLDPC codes is accomplished on the 
Tanner graph through the belief propagation (BP) algorithm 
II20I . II2TI . The BP algorithm exchanges probability vectors of 
length 2™ between variable nodes and check nodes at each 
iteration round £. The BP algorithm mainly consists of 4 parts 
described as follows. 

initialization : We set the iteration round i = 0. Each 
variable node sends the initial message pic""'' = € , 
to each adjacent check node c where v ~ 1,2,...,N and 
c ^ 1,2,..., A/. For general NBLDPC codes (T = 1) 
and multiplicatively repeated NBLDPC codes (T > 2), the 
probability p[?^ (x) is initialized as follows 

p^°\x) =Pr(X„ =r„x|n 

T 

^ Pr(X(tjv+t, = nN+vX\YtN+v = ytN+v), 

t=l 

where ^ represents the normalized constant so that 
'^^i=QP(^2'n-jPv {x) = 1. Note that we defined r^, = 1 e 
GF(2™) for V = 1,...,N. 

Let Vp denotes the set of punctured variable nodes. For 
punctured NBLDPC codes and v ^ 1,2, . . . , N, the probabil- 
ity pi,"' (x) is initialized as follows 



Pi"\x) 



2^ I' e .3. 

Pr(X^ = x\Yy = Hv) otherwise. 



check to variable : For each check node c = 1,2, . . . , M, let 
dc be the set of adjacent variable nodes of c. The check node 
c sends the following message G to each adjacent 
variable node v € dc 



Pi'Kx) =pi'](h-}x) for X e GF(2"), 
pii+'\x) = pii+'\hc.x) for X € GF(2'"), 



where pi ® P2 G is a convolution of pi e which 
can be expressed as follows 



{Pl ®P2){x) 



E 

y,2GGF(2'" 
x=y+z 



Pi{y)p2{z) for X e GF(2'"). 



The convolution appeared above can be efficiently calculated 
via FFT and IFFT 121]. Increment the iteration round as £ := 

£ + 1 

variable to check : For each variable node v = 1, 2, . . . , A^, 
let dy be the set of adjacent check nodes of v. The message 
p^vc € GF(2"') sent from v to c is computed as follows 

pi'}{x) = epi°^(a;)nc'eOA{c}Pc'!,(2^) for ^ G GF(2™), 
where ^ is the normalized constant so that 

tentative decision : For v = 1,2, . . . , N, the tentative decision 
of the w-th symbol is given by 



Letx= 



arg max p^'^ {x) TT p^^' (x) , 

xeGF(2") ^gg^ 

, ij^' ) be the estimated codeword of 
iteration round I. The decoder stops when the maximum 
iteration ^max is reached or Hx = e GF(2™)^^. Otherwise 
repeat the latter 3 decoding steps. After the decoder stops, R 
sends x as the codeword in MAC mode. 

in. Relay Channel 

This paper deals with the time-division half-duplex re- 
lay channel which employs the decode-and-forward protocol. 
Throughout this paper, we refer to the time-division half- 
duplex relay channel with the decode-and-forward protocol 
as the relay channel for the sake of simplicity. 

A. Time-Division Half Duplex Relay Channel 

In half-duplex relaying, a relay cannot transmit and receive 
simultaneously in one time slot so the concept of time-division 
is introduced |[T3l . ll30l . In our model, we assume that the 
relay (R) places on the direct line between a source (S) and 
a destination (D) as shown in Fig. |4] The distance between 
S and D is normalized to 1 and let d denote the distance 
between S and R. The transmission in relay channel takes 
place over two time slots of normalized duration. In BC mode, 
S sends information to R and D within the first time slot of 
duration t. We refer a parameter t as time-sharing factor. In 
the second time slot of duration = 1 — R and S send 
information to D and we refer to the operations in this time 
slot as MAC mode. Figure |5] illustrates the transmission in 
the time-division half-duplex relay channel. Parameters h^D, 
hsK and /ird generally represent the channel effects such as 
fading, shadowing and path loss between two terminals. In this 
paper, we consider only the large scale path loss. Under this 
circumstance, |/isDp = 1, |/isRp = ^ and \hmf 
where a is the path loss exponent. 

For a fair comparison with direct transmission (the commu- 
nication that consists of S and D), the transmitted powers of S 
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Fig. 4. A relay system with a source, a relay and a destination denoted by 
S, R and D respectively. The relay is on the direct line between S and D. 
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Fig. 5. The transmission in the time-division half-duplex relay channel. 

and R must obey an average global power constraint given as 

El 

e : iPs,BC + t'iPsMAC + Prmac) < (6) 

where Ps,bc is the average transmitted power of S in BC 
mode, Ps.MAC ^nd Pr^mac ^6 the average transmitted powers 
of S and R in MAC mode respectively. We define P as the 
total transmission power of the relay system. The alphabet in 
subscript is used to describe the position of variables under 
consideration. 

(i) 

For i-th transmission, let Sg g(-, be the modulated binary 

signals transmitted from S in BC mode, let Sg'j^^^-, be the 
modulated binary signals transmitted from S in MAC mode 

(i) 

and let mac '■^^ modulated binary signals transmitted 
from R in MAC mode. In the similar way, yl)%Q, Vdmac 
Vr BC defined as the received signals at given positions 
and modes. From the notations above, the received signals at 
R and D in BC mode can be written as 

2^D,BC = ^SDSg + 'T-D,BC' 
2^R,BC = ^SRSg + BC 

In MAC mode, the received signal at D is given as 

WD, MAC ~ "SD^s.MAC ^ "■M''R,MAC ^ "d.MAC '^^> 

where u^^'-qq, '^r mac "-d^mac AWGN with zero mean 
and unit variance at given positions and modes. In this paper, 
the relationship between the total transmission power and the 
signal to noise ratio (SNR) of the relay system is defined as 

follows ini. 

SNR =iPs,BC+i'(Ps,MAC + PR,MAc) ,.r.. 

= 2i?,.Eb/No, ^ ' 

where i?r is the coding rate for the relay channel and Eb/No 
is the the energy per bit to noise power spectral density ratio. 
Since the power of noise is normalized to unity, the SNR at 
S in both BC and MAC modes can be expressed as follows 

SNRd3c = |/isd|'Ps,bc, (11) 

SNRd^mac = ( hsu\/Ps,MAC + hw\/ Prmac ) • (12) 



B. Decode-and-Forward Protocol 

We summarize the common parameters before describing 
the decode-and-forward protocol. Let K denote the informa- 
tion length in symbol. Let N denote the length of the code used 
at S for BC mode. Let N' denote the length of the codes used 
for MAC mode. The relationship between the time-sharing 
factor and the lengths of codes used in BC and MAC modes 
can be expressed as -§7 ~ jj- The overall coding rate of the 
relay channel is given by = pf^j^, 

In BC mode, S firstly encodes the K information symbols 
with a code of length N and of rate i?s,BC = The N 
coded symbols are transmitted to both R and D. R decodes 
its received signals and the estimation obtained from decoder 
will be further used in MAC mode. Meanwhile, D just stores 
the received signals. In MAC mode, R transmits another set 
of N' symbols to D. At the same time, S also simultaneously 
sends N' symbols to D. Finally, D can decode the received 
information transmitted in BC mode at lower coding rate 
^D.MAC = N+N' ^ since it receives the additional 
information transmitted in the MAC mode. This implies that 
the computational complexity of decoder at D is much higher 
than that of decoder at R since the number of parity symbols 
at D is more than the number of parity symbols at R. This is 
more crucial when the overall relay system operates at low 
coding rate. In this paper, we will develop coded relay system 
which R and D can decode their received signals with almost 
the same computational complexity. 

C. Achievable Rate 

For a Gaussian relay channel, a general time-division half- 
duplex relay channel with decode-and-forward protocol can 
achieve the following rate ifTTl . 

n= sup mm {tC{w) + t'C{x),tC{y)+t'C{z)}, 
e,o<t,r<i 

(13) 

where 

W =\hsR\^ Ps.BC, 

X =(l - r^) I/isdI^ Ps.mac, 
y =|/isd|^Ps,bc, 

z =|^sDr Ps.MAC + I^rdI^ Pr.mac 

+2r\J\hsTi\^ |/ird|^ Ps,macPr,mac, 
C(SNR) = log2(l + SNR), 

and Q is defined in (|6]l. The parameter r represents the 
correlation between Ss.mac and sr.mac in MAC mode. The 
function C(-) is known as the Shannon formula used to cal- 
culate capacity of Gaussian link. In order to achieve TZ given 
in ( fT3] l, a joint optimization is needed over the correlation 
in MAC mode, time-sharing factor and power allocation. 
Figure |6] compares the capacity of the direct transmission 
with the achievable rate of a relay channel with d = 0.5 and 
Q = 2. In this figure, we plotted the coding rate versus Eb/No 
which is the normalized SNR. Details of these calculations 
are fully presented in iflTl . lISTl . This figure clearly shows 
that the achievable rate of the relay channel is much higher 
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than the capacity of direct transmission especially for low 
SNR. Therefore, we focus our attention on designing low-rate 
NBLDPC codes for the relay channel. 
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■ Relay : Gaussian input 

■ Relay : BPSK input 
Direct : Gaussian input 
Direct : BPSK input 




Fig. 6. The achievable rate of a relay channel for d = 0.5 and a = 2 obtained 
from optimal relay parameters and the capacity of direct transmission. 



IV. Low-Rate Relay Coding Strategy 

In this section, we propose a novel low-rate NBLDPC 
coding strategy for the relay channel. The proposed coding 
strategy is analogous to constructing Ct of the multiplicatively 
repeated NBLDPC code. For this strategy, the relationship 
between time sharing factor t and repetition parameter T 
are t = 4. and hence t' = 1 - ^. For a;,„ e GF(2'"), 



^ and hence t' 
w = 1, 2, . . . , iV and i 1, 2, . . . , m, R, let ^a^g be 
the random variable with the corresponding received signal 



((^^— l)m+i) 



^R.BC 

and K 



We also define the random variables X- 



((i.-l)m+i) 



((v — l)m+i) 
D,MAC 



D,BC 



and these variables are interpreted similarly 
to the definition of ^Jj^^c 

using the decode-and- 
forward protocol, we describe the low-rate NBLDPC coding 
strategy for the relay channel as follows. 

1) Encoding in BC mode: By using an NBLDPC code of 
rate i?s,BC ~ ^ C^i, S encodes K information symbols 
into codeword x = {xi, X2, ■ ■ ■ , x^) G GF(2'")^. S sends x 
to both R and D. For this strategy, S also stores the codeword 
X € Ci for the next transmission in MAC mode. Subsequently 
X is mapped to a binary sequence of length mN bits. Then 
the BPSK modulated signals of the binary sequence are sent 
to R and D at the same time. 

2) Decoding in BC mode: R decodes its received signals 
by using the BP algorithm (described in section II-D). R first 
calculates p'^\x) from J/r^bc^^™''^^^' 



■ ' Vrbc follows. 



i=l 
m 



where (xi, . . . , a;,„) e GF(2)™ is the binary representation 
of a; G GF(2'"), and is the modulates binary signal 
corresponding to Xi. Then the decoder at R is initialized 
pi^"* (x) as the input of the BP algorithm. 

After BP decoding, R produces an estimated codeword Xr G 
GF(2™)^. Meanwhile, D only stores its received signals for 
the future decoding at the end of MAC mode. 

3) Encoding in MAC mode: The encoding processes for 
both S and R in this mode are analogous to constructing 
the code Ct from multiplicative repetition. By using x = 
{xi,X2, ■ ■ ■ ,xn) G Ci stored in BC mode, S produces a 
codeword of code Ct by multiplicative repetition as described 
in Section II-B. 

S then sends only the multiplicative symbol xs = 
{xn+1, xn+2, ■ ■ ■ ,xtn) of length N' = {T—1)N symbols to 
D, where X(^t-i)N+v r^t-i)N+vXv, for w = l,...,N,t = 
2, . . . , T, and ^at+i, . . . ,rTN are randomly chosen from 
GF(2™) \ {0}. 

In the same way, as S produced xs from x, R produces 
xr e GF(2™)('^-i)^ from Ar e GF(2™)^. We note that S 
and R send the same coded symbols Xg = Xr if R, in BC mode, 
did successfully decode the codeword, i.e., x = x. This means 
the additional parity symbols sent from S and R in this mode 
are fully correlated (r = 1). Thus, the additional (T — 1)A^ 
symbols are sent to D. 

4) Decoding in MAC mode: At this step, D has 
two received signals which are j/^^g^^, . . . , yp"|3c^ 

^D^MAC^^ • ■ • ' 2/d'mac' which correspond to the overall TN 
coded symbols or equivalently mTN coded bits. Therefore, 
the overall coding rate of the relay channel is i?r = = 

^^^Y^- The received signals i/d^bc • ■ • ' v'obc correspond to 
the coded symbol of length N transmitted in BC mode, 
whereas the received signals ^^"mac^^ ■ • ■ ' ^d'mac^ corre- 
spond to the coded symbol of length [T ~1)N transmitted in 
MAC mode. Let X,, be the random variables with realization 
x^ where represents the overall coded symbol transmitted 
from both BC and MAC modes and v = 1,2,..., TN. Let 
Yy be the random variables with realization ?/„ which is the 
received value from both BC and MAC modes. The probability 
of transmitted symbol Pr(X„) is assumed to be uniform. From 



7/^™^^ and 



first calculates pl"^ {x) for v 



1, 



' ^DrMAC^ *e decoder at D 
. , iV as follows. 



((ij — l)m+z) 



D,BC 

for i 



{(v-l)r. 

Vd.bc 
. ,m) 



T-1 



Y[ FriXtN+v = TtN+vX I 



({tN +v)m+i) 



T-1 



1 rn 

nn 

t=i i=i 

where st 



D,MAC 
1 



({tN+v)m+i) 
tfD.MAC 



for i ~ 1, . . . ,m) 



'"^P\ 9 ItfD.MAC 



is the modulated binary signal correspond- 
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ing to the i-ih bit of the binary representation of 
rtN+vXv e GF(2™), and ^ is the normalized constant so that 
Sa;6GF(2'") ~ ^- Then the decoder at D is initiaUzed 

p^\xv) as the input of the BP algorithm. After BP decoding, 
D produces an estimated codeword xd E GF(2'")^. 

The initial calculation of p^'' (x) can be viewed as being cal- 
culated by combining two received signals j/n'g(j and yn^MAC- 
After the two received signals are combined, the decoder 
performs the BP algorithm on the Tanner graph of Ci with N 
variable nodes. 

In summary, the proposed low-rate NBLDPC coding strat- 
egy for the relay channel is depicted in Fig. |7] We claim 
that the proposed coding strategy is very simple due to the 
following two reasons. 1) For encoding in MAC mode, the 
encoder at R uses only multiplicative repetition which is much 
simpler than encoders of BLDPC codes. The multiplicative 
repetition requires only {T—1)N multiplications over GF(2™) 
for encoding. 2) The Tanner graphs of decoders at both R and 
D are the same, i.e., the Tanner graph of Ci. The difference 
of decoding at R and D is only the initialization. And the 
computational complexity of the decoders are almost the same 
even if the coding rate at D was lower than that of R. 



Source 



information of lengtii K symbols 



xeC, 



xeC, 



Relay 



Destination 



MR 



xgC, 



i/R,i 



MR 



i/D,BC 



i/D.MAC 



combine 



• BC mode 

• MAC mode 



Decoder 



▼ 



Fig. 7. Block diagram of the proposed low-rate NBLDPC coding strategy 
for the relay channel. In this figure, the block labeled with "MR" represents 
multiplicative repetition as described in Section IV. The block labeled with 
"combine" represents the mathematical operations according to the initializa- 
tion step of decoding at D. The Tanner graphs at the decoders at R and D are 
the same. 



V. Numerical Results 

In this section, we demonstrate the decoding performance 
of NBLDPC coded relay system. For all results, the NBLDPC 
codes defined over GF(2^) with regular degree profile {d^ = 
2,dc) are employed. We choose 7ti = 8 for its good per- 
formance and also the computer-friendliness of byte. The 
FFT-based BP algorithm i211 with the maximum iteration 
^max = 500 is used for decoding. The BPSK modulation 
was used for all the links in the relay channel. We set the 
attenuation exponent to a = 2 and the source-relay distance 
to d = 0.5. We denote the bit error rate and frame error rate 
as BER and FER respectively. 

First, we present the BER performances of the proposed 
low-rate coding strategy for the relay channel. The power al- 
location for S and R are static and fixed as follows: Ps.bc = "j, 



n 



S.MAC 



and Pr. 



MAC 



A (dv — 2,dc = 3)-regular 



NBLDPC codes of rate i? = | over GF(2^) are chosen as the 
mother code Ci . We choose the repetition parameter T = 2 for 
encoding in MAC mode. This corresponds to the equal time 
sharing constraint since t = y = ^. By using this strategy, 
the overall coding rate is given by = ~ ~ \- 

Figure [8] shows the BER performances of the proposed 
three low-rate NBLDPC codes for the relay channel. The three 
NBLDPC codes have information length k G {56, 192, 1024} 
bits, respectively. The threshold value for the proposed coding 
strategy with i?,- = ^ is calculated by the Monte Carlo density 
evolution. The Monte Carlo density evolution originally devel- 
oped in Jill is a method to find the threshold of NBLDPC 
codes in the limit of very large codeword length. It can be seen 
that the proposed codes theoretically perform very close to 
the achievable rate. The gap of threshold from the achievable 
rate is only 0.3 dB. We can see from the figure that BER 
curve improves by increasing of information length (also the 
codeword length). At BER of 10^^ and = i, it can be 
seen from the figure that the proposed NBLDPC code of 
information length k = 1024 bits performs within 1.2 dB 
from the achievable rate. Moreover, the coding gain is about 
1.1 dB beyond the capacity of direct transmission. 

We explain how the proposed strategy can be efficiently 
used for the rate-compatible relay system by investigating the 
performances of relay systems with i?, = i obtained from 
puncturing the codes used in relay systems with i? = i. 
By using the recoverable step discussed in Section II-C, 
we puncture the (2,3)-regular mother NBLDPC code of rate 
i? = i to obtain an NBLDPC code of rate i?s,BC = 5- 
For example, we obtain a (A^ = 48, K = 24) NBLDPC 
code by puncturing {N = 72, K = 24) NBLDPC code. 
This means that we cancel Np = 24 parity symbols of 
(TV = 72, K = 24) NBLDPC code. With repetition parameter 



T = 2 at R, the overall coding rate is R^ = = i. The 

BER performances of the proposed strategy with punctured 
NBLDPC codes (dashed curves) are shown in Fig. [8] Although 
the codes are punctured but the NBLDPC coded relay system 
of rate R = j still obtain a coding gain around 0.75 dB beyond 
the capacity of direct transmission. 

As mentioned in ifTTl , the optimized BLDPC codes need 
the concatenation with outer BCH or RS codes to guarantee 
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-3 -2 -I 



E^/N„(dB) 

Fig. 8. The BER performances of NBLDPC coded relay systems for low- 
rate case. The parameters in parentheses represent the overall transmitted 
bits and the information bits respectively. In this figure, the word "Relay" 
represents achievable rate of relay channel and the word "Direct" represents 
the Shannon's capacity of direct transmission. The threshold value for i?r = ^ 
is represented by "Threshold". The black curves labelled with "puncturing" 
represent the rate-compatible relay systems with /?r = j obtained by 
puncturing the codes of rate iJr = 4 . 



the good FER performance. The advantage of using {2,dc)- 
regular NBLDPC codes over the direct transmission is that 
the regular NBLDPC codes with dv = 2 have a good FER 
performance at both waterfall and error region 1331 . 1341 . 
Figure |9] shows the FER performances of NBLDPC codes 
for the relay channel. The proposed low-rate codes exhibit 
excellent FER performances at both waterfall and error region 
without concatenating BCH or RS codes. The punctured 
NBLDPC coded relay systems with R = j (dashed curve) 
exhibit very good FER performances. 




E^/N„(dB) 

Fig. 9. The FER performance of NBLDPC coded relay systems for low-rate 
case. 



VI. Conclusion 

We propose in this paper the low-rate NBLDPC coding 
strategy for the decode-and-forward half-duplex relay chan- 
nels. The advantages of the proposed strategy can be listed 
as follows; 1) The waterfall performance within 1.5 dB from 
the achievable rate can be obtained from the moderate length 
NBLDPC codes. The significant coding gain is also achieved 
over the Shannon's limit of the direct transmission. 2) We do 
not need the optimization process for optimizing the degree 
profile of LDPC codes. We just employ the (2, dc) -regular 
NBLDPC defined over GF(2*^) as the channel code. 3) The 
FER performances obtained from the (2, dc)-regular NBLDPC 
codes are excellent even the outer BCH or RS codes are 
not employed to lower the error floors. 4) The relay and the 
destination can decode with almost the same computational 
complexity since both the relay and the destination use the 
same Tanner graph for decoding. 5) The encoding processes 
for both the source and the relay in MAC mode are very simple 
since we use the multiplicative repetition which requires only 
(T — l)iV multiplications over GF(2'"). 6) The proposed 
strategy is applicable for rate-compatible relay systems. By 
adapting the repetition parameter T at the relay or puncturing 
the codes at the source, easy HARQ is possible for very noisy 
relay channels. 
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